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ABSTRACT 

Aims. We present and analyse late-time observations of the type-Ib supernova with possible pre-supernova progenitor detection, 
iPTF13bvn, taken at ~300 days after the explosion, and discuss these in the context of constraints on the supernova’s progenitor. 
Previous studies have proposed two possible natures for the progenitor of the supernova, i.e. a massive Wolf-Rayet star or a lower- 
mass star in close binary system. 

Methods. Our observations show that the supernova has entered the nebular phase, with the spectrum dominated by Mg I]/I/i4571, 
[O I]/l/16300, 6364, and [Ca 11]AA7291, 7324 emission lines. We measured the emission line fluxes to estimate the core oxygen mass 
and compare the [O I]/[Ca II] line ratio with other supernovae. 

Results. The core oxygen mass of the supernova progenitor was estimated to be <0.7 Mq, which implies initial progenitor mass not 
exceeding ~-15 - 17 Mq. Since the derived mass is too small for a single star to become a Wolf-Rayet star, this result lends more 
support to the binary nature of the progenitor star of iPTF13bvn. The comparison of [O I]/[Ca II] line ratio with other supemovae also 
shows that iPTF13bvn appears to be in close association with the lower-mass progenitors of stripped-envelope and type-II supemovae. 

Key words, supemovae: general - supemovae: individual: iPTF13bvn - stars: massive 


1. Introduction 

Supemovae (SNe) type Ib and Ic, wbicb do not show hy¬ 
drogen lines in tbeir spectra, are thought to be produced 
by the core-collapse of massive stars which have lost their 
outer hydrogen envelope. Hydrogen-deficient Wolf-Rayet (WR) 
stars have been proposed a s the p ossible progenitors f o r thes e 
SNe (iMaeder & Leaueu^ Il982l: lEnsman & Woosle^ Il988h . 
However, thus far there have been no strong c onfirmation of 
WR st ar as SN Ib/c progenitor. As summarized bv lEldridee et al.l 
(l2013l) . examinations of Hubble Space Telescope (HST) pre-SN 
archival images of SNe Ib/c in nearby galaxies has not suc- 
ceeded in detecting the progenitor stars. lYoon et ^ (l2012h ar¬ 
gued that type-Ib/c SN progenitors are faint in optical bands dur¬ 
ing the pre-SN stage thus direct detections have not been provid¬ 
ing strong constraints on the progenitor luminosity and mass. 


* Based on observations obtained at the Southern Astrophysical 
Research (SOAR) telescope, which is a joint project of the Ministerio 
da Ciencia, Tecnologia, e Inovajao (MCTI) da Republica Federativa do 
Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the 
University of North Carolina at Chapel Hill (UNC), and Michigan State 
University (MSU); Chilean Telescope Time Allocation Committee pro¬ 
posal CN2014A-9L 
** e-mail: hanin@das.uchile.cl 


ICao et all d2013ll reported the discovery and early evolu¬ 
tion of iPTE13bvn, a type-Ib SN discovered by the intermedi¬ 
ate Palomar Tran sient Factory (iPTF) survey (iRau et al.L 120091: 
iLaw et aQ l2009h in the host galaxy NGC 5806. It was found 
in the HST archives that the SN position coincides with a blue 
source, within 2cr (80 mas, « 8.7 pc) positional uncertainty of 
the SN. The photometry of the object was shown to be consis¬ 
tent with a WR star, which w as suggested as the progenitor of 
iPTF13bvn. lGroh et al.l(l2013h proposed that their model of WR 
star with initial mass of 31 - 35 Mq could reproduce the ob¬ 
served photometry of the object, the candidate progenitor star 
of iPTF13bvn. If the disappearance of the candidate progenitor 
star is confirmed by future observations, this would be the first 
detection of a progenitor of a type-Ib SN. 

An alternative scenario was proposed by iBersten et al.l 
(12014 . in which the progenitor star was modeled to be in a 
close binary system. Their hydrodynamic model shows that the 
primary star in this system evolved and lost its hydrogen layer 
due to binary interaction, and exploded as a ^3.5 Mq pre-SN 
star, corresponding to the progenitor initial mass of 20 Mq. This 
model successfully reproduced the observed light curve, absence 
of hydrogen, and photometry of the pre-SN candidate progenitor 
star. The notion that the progenitor star of iPTF13bvn could not 
have been a massive WR star, but instead a lower-mass star in bi- 
nary system, has also been suggested by iFremling et alJ d20l4 
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Table 1; Observing log. 


Object 

Exp. time (s) 

Airmass 

Seeing (”) 

iPTF13bvn spec. 

5100 

1.4- 1.2 

1.2-0.9 

HR5501 spec. 

1 

1.2 

0.9 

iPTF13bvn BVRr 

300 

1.2- 1.3 

p 

1 

iPTF13bvn U 

1800 

1.3 

1.2 

PG1528 BVRr 

1-3 

1.4 

1.3 

PG1528 U 

30 

1.4 

1.3 

PG1525 BVRr 

2-5 

1.3 

1.2 

PG1525 U 

30 

1.3 

1.2 


Note: HR5501 is a spectrophotometric standard star; PG1525 and 
PG1528 are Landolt photometric standard fields. 


and lSrivastav et al.l (12014 K Eldridge et all (1201 5h argued that the 
photometry of icao et al. ( 2013h underestimated the brightness of 
the candidate progenitor object by up to 0.7 mag, which brought 
them into the conclusion that the observed spectral energy distri¬ 
bution does not match a typical single WR star, but rather more 
consistent with their binary models. 

In this paper we present optical photometry and spectroscopy 
of iPTF13bvn obtained during its nebular phase, at 306 days af¬ 
ter the explosion date, from which we derive an estimate of the 
initial mass of the SN progenitor star. 


2. Observations and data reduction 

The observations of iPTF13bvn were carried out using the 4.1 
m Southern Astrophysical Research (SOAR) telescope stationed 
at Cerro Pachon, Chile at an altitude of 2700 m above sea 
level, starting in the night of 18 April 2014 (Chilean local time) 
in visitor mode. This observation date corresponds to ~306 
days after the assume d explosion time of JD = 2456459.24+0.9 
dBersten et all l2014ll . or +290 days past the time of maximum 
light. 

We used the Go odman High Throughput Spectrograph 
dClemens et al.L l2004ll for both photometry and spectroscopy. 
The Goodman HTS employs a Fairchild 4kx4k CCD with a 
physical scale of 0.15”/pixel. The unvignetted field of view size 
of the instrument is around 7.2’ in diameter. We used the 2x2 
binning mode for both photometry and spectroscopy. During the 
night, the sky was clear but not photometric, with variable see¬ 
ing. Table [U shows the observation log for the night. The obser¬ 
vation of iPTF13bvn was started with spectroscopy, followed by 
photometry. Sp ectrophotometri c standard star and photometric 
standard fields dFandold Il992h were also observed during the 
night. 

We used the 300 lines/mm grating with GG385 blocking fil¬ 
ter and 1.03” slit mask, positioned at the parallactic angle, for 
our spectroscopic observation of iPTF13bvn. With this instru¬ 
ment configuration, the dispersion for 2x2 pixel binning is 2.6 
A/pixel with wavelength coverage 3600 - 8800 A. We took a 
1500 s exposure and 2x1800 s exposures of iPTF13bvn, result¬ 
ing in a total of 5100 s exposure time for the object. We used 
the 3” slit for the spectrophotometric standard star. An internal 
HgAr arc lamp was used for wavelength calibration, and inter¬ 
nal quartz lamp was used for spectral flat calibration. From the 
arc spectra we measured the effective full width at half maxi¬ 
mum (FWHM) resolution of the spectra to be xH .6 A around the 
wavelength 7000 A. The raw data was reduced using IraiQ in a 


standard slit-spectroscopy reduction manner, with the twodspec 
and ONEDSPEC packages. Wavelength and flux calibrations were 
also performed. 

Photometry was carried out in UBVRc bands. Figure[T]shows 
the position of the SN within its host galaxy. The SN is clearly 
detected in the Rc band, as well as in B and V although not 
as clearly. In the t/-band, the SN is not detected within 2cr of 
the background noise. The integration time was 300 s for the 
BVRc filters, and 1800 s for the U filter. Point-spread function 
(PSF) photometry was performed to the BVRc images and aper¬ 
ture photometry was done to the t/-band image using aperture 
value equal to the seeing size, after standard image reduction 
procedures with I raf. PSF photo metry was performed using the 
DAOPHOT package (IStetsonLll987h . and aperture photometry using 

APPHOT. 


3. Results and discussions 

3.1. Spectroscopy 

Our observation was conducted at 306 days post SN explo¬ 
sion. At this age, the SN is expected to have already en¬ 
tered the nebular phase, which is confirmed by the observed 
spectrum (Figure |2]l. The continuum of the spectrum is very 
weak, with strong emission lines of [O I]/l/l6300, 6364 and [Ca 
II]/l/l7291, 7324 dominating. In addition to the [O I] and [Ca 
II] lines, the spectrum also exhibits a noticeable Mg I]/14571 
semi-forbidden intercombination line. These spectral features 
are commonly shown by stripped-envelope core-collapse SNe 
at ne bular stage (see Figure [Hi. Nebular spectra of s even other 
SNe (iModiaz et al.Ll201^lM^da et alilT0()7bll20d^ are shown 
alongside iPTF13bvn for comparison in Figure|2l 

The nebular spectrum of iPTF13bvn looks remarkably sim¬ 
ilar to the type-IIb SN 2008ax at +307 days post-maximum, 
despite the noticeably weaker signal-to-noise ratio. SNe lib are 
thought to be closely related to SNe Ib - both subclasses show 
prominent helium lines, but the hydrogen lines present in type- 
IIb SNe are observed only during the early phase. Thus, it is 
possible that many SNe lib might have been misidentified as a 
type-Ib. Indeed it has been sugg ested that the Ilb or Ib cla ssi- 
fication may be time-dependent (iMilisavlievic et all l2013h . In 
the case of iPTF13bvn, hydrogen was not detected in spectra as 
early as +2 . 3 day s after the explosion, hence the Ib classification 
(ICao et al.L l2013l) : howeyer, see the discussions about possible 
Ha emission in subsection l3.1.3l 

3.1.1. The [O I] line profile 

The [O I]/l/l6300, 6364 emission line of iPTF13byn eyidently 
does not show a single-peak, Gaussian-like profile. It exhibits a 
rather boxy, flat-topped profile, that to some extent is sugges- 
tiye of a double-peaked emission. In Figure [3] we compare the 
line profile of the [O I]/l/l6300, 6364 emission in iPTF13byn 
to that of SN 2008ax at the same phase. The spectrum of SN 
2008ax is the same as the one presented in Figure|3 Both spec¬ 
tra show [O I] FWHM of -5000 km s '. While the [Ca II] line 
of iPTF13byn also show a similar yelocity FWHM of -5000 
km s“', in the case of SN 2008ax the FWHM is lower, -4000 
km s '. These are within the likely c ore yelocity yalues of 
- 100 0 - 5000 km s ' as suggested by iFransson & Cheyaliej 

dHH. 


' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 


for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 1: Sections of UBVRc images of iPTF13bvn, with the SN indicated. 


Mg I] 4571 [0 l] 6300,6364 [Co ll] 7291,7324 



4000 5000 6000 7000 8000 9000 

Rest wavelength (A) 

Fig. 2; Spectrum of iPTF13bvn at 306 days past explosion (+290 days after maxim um light), com pared to seyeral strip ped-enyelope 
core-collapse SNe at similar ages from iModiaz et al.l (l2014ll . lMaeda et ahl (l2007bi SN 2005bf), and iMaeda et al.l (1200^ SN 2004gq). 
The phase of the spectra in days are with respect to the time of maximum light. The prominent emission lines originated from the 
SNe are indicated with yertical grey lines. The spectra are normalized to the peak of the [O I] line of iPTF13byn and shifted in flux 
for presentation purposes. The narrow emission lines originate from the underlying H II regions. 
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Fig. 3; Line profile of [O III]/1/16300, 6364 of iPTFBbvn com¬ 
pared to SN 2008ax. The dashed vertical line corresponds to zero 
velocity with respect to 6300 A. The flux of the spectrum of SN 
2008ax has been normalized and shifted. 


In Figure ID we present the profiles of the [O I], [Ca II], and 
Mg I] emission lines. Here we plot the observed line profile 
of Mg I] with full line, and additionally a synthetic line pro¬ 
file with dash-dot line. This synthetic prohle was generated to 
simulate a doublet feature for Mg I], to allow direct compar- 
is on with the [O I] doublet . We follow the method described 
in iTaubenberger et ahl (|2009|) . in which the Mg I] line is scaled 
down to 1/3 of its initial intensity then shifted by -1-46 A and 
adde d to the original spectrum to generate the synthetic dou¬ 
blet. ITaubenberger et aH (l2009l) found that for most objects in 
their sample (~ 65%) the [O I] and Mg I] lines show remarkably 
similar profiles. It is interesting that even with a synthetic dou¬ 
blet, the red wing of the Mg I] line in iPTF13bvn seems to be 
depressed in comparison with the profiles of [O I] and [Ca II] 
although there is the caveat of low signal-to-noise ratio. The 
strength of this line is also considerably stronger compared to 
some other SNe as apparent in Figure |3 The Mg I]/[0 I] flux 
ratio in iPTF13bvn is around 0.85, whereas for S Ne 2002ap, 
2008a x, and 2007Y the value is around 0.2 - 0.3. iFolev et S 
(l2003ll noticed the exceptional strength of this Mg I] line in SN 
2002ap at late times and that it is growing with time compared 
to the strength of [O I] and [Ca II] lines. As time passes we are 
seeing deeper into the progenitor core and reach the Mg-O layer, 
thus it is expected that this line grows stronger at later times. For 
the [Ca II] line in iPTF13bvn, the line profile appears to be less 
asymmetric compared to the other two lines. These might indi¬ 
cate different geometry and distribution of the emitting materials 
for[0 I], [Call], and [Mg I]. 

The profile shape of the [O I] line in iPTF13bvn suggests 
that there is some degree of asphericity in the ejecta. While it 
is not very clearly double-peaked, the line profile is consider¬ 
ably d ifferent to the more Gaussian-like [Ca II] line. lMaeda et al.l 
(I2OO8I) showed that double-peaked [O I] 4/16300, 6364 emission 
commonly observed in late-time spectra of core-collap s e SNe 
signifies asphericity in the SN explosion. iTanaka et al.l (l2009t) 
also reported that the [O I] line in the nebular spectrum of the 
type-Ib SN 2008D can be explained by a bipolar explosion wit h 
a torus-like distribution of oxygen. ITaubenberger et aP (1201 ll) . 
who also observed the similar feature in SN 2008ax, argued 


•< 

I 

E 

u 

cr 

(D 

o 



- 6000-4000 -2000 0 2000 4000 6000 

Velocity (km s ) 


Fig. 4; Line profiles of [O I], [Call], and Mg I], in velocity space. 
Zero velocities are with respect to 6300 A, 7308 A, and 4571 A, 
respectively. The fluxes of [Ca II] and Mg I] lines are normalized 
to the [O I] line. The synthetic Mg I] line profile is represented 
with dash-dot line (see text for details). 


that a spherical oxygen distribution with a significant clump or 
torus, or even an aspherical distribution of ^®Ni that excites the 
oxygen, may explain the observed feature. On the other hand, 
the asymmetry in the line profile could be attributed to ab¬ 
sorption in the interior, which may be caused by dust, or high 
opa city associated with clumping or high density mater ial (see 
e.g. ITaubenberger et al.Ll200^lMilisavhevic et aill201(ll) . In this 
case, the observed double-peaked profile may also not neces¬ 
sarily reflect a torus or elongated-shell geometry of the oxygen 
emitting region, but simply the dou blet nature of the line as sug¬ 
gested bv iMilisavlievic et alJ (l2010l) . Although it cannot be ruled 
out in all cases, the doublet explanation cannot account for the 
observed ratios of the components, nor for the cases where the 
profile is singly peaked. 

3.1.2. Oxygen and progenitor mass 

Using the flux of the [O I]/l/16300, 6364 emission line, the 
mass of oxygen that produced the line can be estimated. 
Llerkstrand et al.l (l2012h have demonstrated that the strength of 
the [O I]/l/l6300, 6364 line is sensitive to the variation of the 
init ial mass of the S N progenitor star. 

lUomotol (Il986h provided an equation to calculate the mini¬ 
mum oxygen mass (in M© unit) required to produce the emis¬ 
sion, as follows; 

Morgen = 10V([O I])Z)2 exp(2.28/7’4). (1) 

Therefore, using this equation one can calculate the oxy¬ 
gen mass from the measured flux of the [O I]/l/16300, 6364 line 
in erg s“* cm“^, /([O I]), assuming that the distance in Mpc 
D and temperature of the oxygen-emitting region T^, in 10'* K 
units, are known. The distance to the host galaxy of iPTF13bvn 
NGC 5806 derived with Tully-Fisher method is known to be ~26 
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Mpc. The Extragalactic Distance Databas^l (iTullv et all l2009h 
gives the distance to NGC 5806 as 26.2 +1.5 Mpc, and this 
has been revised with a newer value of 2 6.79 ± 0.2 Mpc in the 
Cosmicflows-2 catalog (iTullv et al.Ll2013h . We adopt this newer 
value as th e distance to iPTF13b vn. For reference, ICao et all 
(1201 3h and ISriva stav et ah ( 2014ll adopt the distance of 22.5 
Mpc, while lBersten et al.l (l2014ll uses 25.5 Mpc. From Equation 
[1] it is clear that adopting smaller distance would decrease the 
estimate of the oxygen mass, therefore the progenitor mass. 


The measurement of the emission line flux of [O I]/l/16300, 
6364 was done using the task splot in Iraf with three different 
methods. Prior to the flux measurements, we corrected the spec¬ 
trum of iPTF13bvn for Milky Way and host galaxy extinctions, 
amou nting E(B - V)mw - 0.0447 mag and E{B - V)hnst - 0.17 
mag (Bersten et al.L l2014ll . assuming a standard ICardelli et al.l 
(11989 ) interstellar reddening law with Ry = 3.1, and removed 
the spectral continuum by fitting a first-order cubic spline func¬ 
tion. Measurement of the line flux by direct integration gives 
3.46 X 10“'^ erg s“' cm“^, by fitting a single gaussian function 
3.96 X 10 erg s * cm“^, and by deblending the profile into two 
gaussians 3.70 x 10“'^ erg s“' cm“^. We use the average of the 
three measurements, 3.71 ± 0.25 x 10“'^ erg s“' cm“^ as the flux 
of the oxygen line, to be used in equation [T] A possibly sign if- 
icant source of uncertainty in this method of lUomot^ (119861) is 
the time evolution of the strength of the [O I] line as well as the 
temperature T 4 . 


One way to estimate T 4 is by measuring the ratio of 
[O I]/15577 line to [O I] 4/16300, 6364. However, in the case 
of iPTF13bvn we could not detect this weak line at 5577 A, 
which is actually an indication of low temperature (lElmhamdii 
l201lh . The near-infrared Ca II line was also not detected in 
the i PTF13bvn spectrum, f urther indicating the low tempera¬ 
ture. [S^lerman|eraI](jT998l) used this method of [O I]/15577 line 
to estimate the temperature to be < 5000 - 4400 K, which in 
turn gives a n oxy gen mass of 0.11 - 0.21 Mq for SN Ib 1996N. 
ISahu et aTI ( 1201 lb . not detecting the [O I]/15577 in the nebular 
spectra o f SN Ib 2009if, used T 4 - 0.4 to derive oxygen mass of 
1.34 Mq. iFlmhamdi et"^ (|2004 . using the upper limit for the 
[O I]/15577 line, derived oxygen mass of 0.7 - 1.35 Mq for SN 
Ib 19901, for temperatures between 3200 - 3500 K at +237 days. 


Assuming the lower limit of the temperature of the oxygen- 
emitting region to be 3200 K as in SN 19901, while adopting 
the line flux 3.71 + 0.25 x 10“'^ erg s“* cm“^ and distance 
26.79 + 0.2 Mpc, we calculated the minimum oxygen mass in 
iPTF13bvn using Equation dTJ to be 0.33 + 0.03 Mq. This oxy¬ 
gen mass value is relatively low compared to SNe 19901 and 
2009jf, but comparable to SN 1996N if the oxygen temperature 
of iPTF13bvn is assumed to be slightly higher. With higher tem¬ 
perature, Equation ([1]) will yield an even lower value of oxygen 
mass. Again here we note the sensitivity of Equation ([1]) to the 
temperature assumption - increasing the temperature from 3200 
K to 5000 K will drop the derived oxygen mass by a factor of 
~13, while on the opposite direction using very cool temperature 
less than 2500 K will increase the derived oxygen mass by one 
order of magnitude. At any rate, a temperature less than ~3000 
K is not plausible as beyond this value the cooling will shift to 
be dominated by far-infrared fine-structure rather than optical 
lines, resulting in a thermal instab ility and sudden temperature 
drop (iFransson & Chevalieiill98^ . Thus, it is not likely that the 
minimum oxygen mass is significantly higher than the derived 
value of 0.33 Mq 


It is to be noted, however, that Equation[T]provides the mini¬ 
mum oxygen mass responsible for the line emission, which may 
not be equal to the total oxygen mass. The presence of non- 
optically thin materials and clumping is not accounted for, thus 
the total oxygen mass may actually be higher than the derived 
value of 0.33 Mq. iMaeda et al.l (l2007ah derived a total oxygen 
mass of 2.3 Mq for SN 2006aj, 1.3 Mq of which is within 
V < 7300 km s“' and consistent with the luminosity of the 
[O I]/l/16300, 6364 lines. With a conservative assumption that 
~ 50% of the total oxygen mass in iPTF13bvn is responsible 
for the emission line, then the total oxygen mass should not ex¬ 
ceed ~0.7 Mq. The effects of temperature assumption and non- 
emissive oxygen seem to dominate the uncertainty in the oxygen 
mass derivation, comparing with the uncertainties in distance or 
flux measurement. 

In Figure |5] we compare the derived oxygen mass with nu¬ 
cleosynthesis oxygen yields of massive stars of various i nitial 
masses, from iNor noto et all (Il997h . iRauscher et akl (l2002l) . and 
iLimongj & Chieffir i 2003h . It is apparent that the derived oxy¬ 
gen mass of iPTF13bvn is more consistent with what would 
be produced from the explosion of a massive star with mass 
< 15-17 Mq, even if the total oxygen mass is assumed to 
be twice the derived value from the [O I]/l/16300, 6364 line. 
In an extreme case in which the whole S N ejecta (< 2.3 Mq; 
iBersten et al.L 1201 4t ISrivastav et al.L 1201 4l) were assumed to be 
oxygen, which is not realistic, the initial progenitor mass of 
iPTF13bvn would still be smaller than ~25 Mq. For reference, 
the derived oxygen masses of SNe 199 01 and 20091 f are around 
14 - 36 % of th e total ejecta mass (lElmhamdi et al.L 120041: 
ISahu et akl 1201 ll) . This implies that a single Wolf-Rayet pro¬ 
genitor scenario for iPTF13bvn, whose initial mas s must have 
been > 25 Mq at solar metallicity (ICrowtheil l2007l) . is unlikely. 
Indeed, a lower-mass binar y progenitor scenario for iPTF13bvn 
seem s to be more plausible (Bersten et al.L 2014 Fremling et al.L 
l20l4 ISrivastav et al.L 12014 : Eldridge et al. . 20151) . With lower 
metallicity, the WR mass limit would move to even higher mass 
due to the reduced strength of the metallicity-driven wind. In this 
context, the metallicity of the explosion site of iPTF13bvn was 
measured to be slightly less than solar (subsection 13.1.4b . We 
note that nucleosynthesis yields are possibly subject to uncer¬ 
tainties i ntroduced by neglectin g the effects of turbulence in the 
interior (ISmith & Arnettll2014l) . hence this may affect the esti¬ 
mate of the initial mass of the SN progenitor star. However, this 
effect is still not thoroughly studied and quantified. 

We also measured the line ratio of [O I]/l/l6300, 
6364/[Ca 11]AA7291, 7324 emissions in the nebular spectrum of 
iPTF13bvn, and compared it to other stripped-envelope SNe as 
well as several type-II SNe taken from the SUSPECT0 Online 
Supernova Spectrum Archive. This line ratio is known to be in¬ 
sensitive to the density an d temperature, while increase s with in- 
creasing progenitor m ass (IFransson & Chevalieilll98^ also see 
lElmhamdi et al .120041) . The evolution of the ratio of this line with 
SN age is generally small, as shown in an example of SN 2004gq 
in Figured The ratio of [O I]/l/16300,6364/[Ca II]447291,7324 
changes from around unity at phase +258 days to 1.3 at +381 
days, i.e. only ~30 % change in mo r e than 120 days. It is also 
has been shown by lElmhamdi et al.1 (l2004 see their Figure 3) 
that the ratio is almost constant for SNe aged 280 - 400 days 
after the explosion. Figure |7] shows the measured line ratio for 
iPTF13bvn compared to other SNe. It is clear that the SNe II 
do not exceed line ratio of 0.7, while on the other hand there 
is considerable spread for the type-Ib/Ic SNe. While the evolu- 


^ http://edd.ifa.hawaii.edu/ 


^ http://www.nhn.ou.edu/-suspect/ 
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Fig. 6: Evolution of the nebular spectrum of SN 2004gq around 
the [O I] - [Ca II] spectral region, between +258 and +381 days 
post maximum. Spectra has been scaled to m atch the [O I ] line 
at +258 days. Data from iMaeda et al.l (l2008h . 


Fig. 5; Oxygen yield of the explosion of mas sive stars of various 
main sequence m asses ( Nomoto et al.|| 19971: red. iRauscher et all 
l2002t light blue. iLimongi & Chiefhl orange). Horizontal 
line indicates the estimated oxygen mass for iPTF13bvn assum¬ 
ing T - 3200 K, with uncertainty represented by shaded region. 
Dashed horizontal lines denote the oxygen mass values for dif¬ 
ferent temperature assumptions. The estimated oxygen mass if 
50% of oxygen was unseen is indicated with dash-dot line. 


tion of the line ratio may contribute to the scatter in the dia¬ 
gram, it is considerably insignificant compared to the observed 
spreacQ. This spread may be interpreted as indicative of the 
presence of two different populations of SNe Ib/c (shaded re¬ 
gions in Figure|2]l, i.e. those coming from massive single Wolf- 
Rayet progenitors, and low er-mass progenitors in binary systems 
(I Kuncarayakti et all l2Q13h . In this case, iPTF13bvn appears to 
be more associated to the lower-mass progenitors of type-II SNe, 
further strengthening the notion of its progenitor in a binary sys¬ 
tem. It has been suggested that the majority (>70 %) of massive 
stars undergo binary interaction at some point in their evolutions, 
thus implying that a significant fraction of cor e-collapse SN pro - 
genitors must experience such phase as well dSana et al.Ll2012h . 
The conhrmation that the progenitor of iPTF13bvn indeed be¬ 
longed to a binary system may eventually be provided by the de¬ 
tection of the remaining companion by future observations when 
the SN has faded consi derably, such as the case of SN Ilb 201 Idh 
dFolatelli etalll2014bll . 

3.1.3. Possible presence of Ha? 

As can bee seen in Figure |2l the global appearance of the 
late-tim e spectrum of iPTF13bvn c losely resembles that of SN 
2008ax. iTaubenberger et aP (l201 ih reported the presence of a 


We, however, again caution the reader of the caveat that it is pos¬ 
sible that this pattern observed in Figure |7] may change with time as 
different SNe may exhibit different rates of evolution, i.e. some objects 
may evolve more dramatically compared to the others. Note that this 
late-time behaviour in SNe is still not very well observed and studied. 


broad feature redward of the [O I]/l/l6300, 6364 emission line in 
nebular spectra of SN 2008ax. This feature appeared after ~100 
days after the explosion and remained visible up until at least 
phase 358 days. In other stripped-envelope S Ne this feature has 
also been observed, such as in SN Ilb 1993J (iPatat et allll995h 
and SN Ib 1996N (Sollerman et al| |l998l) . In the nebular spectra 
of SN Ib 2005bf (iMaeda et al.Ll2007bh . a broad feature at -6500 
A has been identihed as Ha (see Fig. |2]). Also in the case of 
SN 2010as, a low-velocity SN I lb showing Ib/c charac teristics, 
this feature was also observed dFolatelli et al.L l2014ah . During 
the transitional phase from photospheric to nebular at around 
+ 100 days past maximum light, the profile is complex, possi¬ 
bly indicating the presence of other lines, and appears to be still 
present well into the nebular phase at +309 days. 

Although it is compelling to attribute this feature to late- 
time Ha emission from the SN, as has been discussed by 
ITaubenberger et ahl d201 ih there could be other possible sources 
of this broad feature such as various other elements includ¬ 


ing Fe, and scattering. In Figure 0 we plot the spect r um of 
iPTF13bvn alongside SN 2008ax from iModiaz et all d20l4 
same data shown in Figure |2| around the Ha region. The graph 
is suggestive of the presence of a broad feature redward of the 
[O I] line in the iPTF13bvn spectrum, very similar to SN 2008ax. 
However, due to the possible ambiguity in the identification as 
discussed above and low signal-to-noise ratio in the spectrum 
we restrain ourselves from hrmly associating t his feature with 
Ha em ission from the SN itself. Furthermore, iJerkstrand et al.l 
d2015l) argued that the thin hydrogen envelope in type-IIb SNe 
should not affect the spectrum after around 150 days. They 
attribute the feature redward of the [O I] line in the nebular 
spectra of SNe 1993J, 2008ax, and 2001dh to emission from 
[N II]TT6548,6583. 

We note that the early-phase spectra of iPTF13bvn bear 
similarities with SNe 2007Y and 2010as, which were shown 
to be long to the class of ’’flat-velocity type-IIb” transitional 
SNe dFolatelli et al.L 12014ah . These objects show SN type-Ib/c 
signatures in their spectra, but also the presence of hydro¬ 
gen, and characteristic flat velocity evolution between 6000 and 
8000 km s“5 The velocity evolution of iPTF13bvn, however, 
does not show this characteristic flat evolution (see Figure 9 of 
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Fig. 7: [O I]/i^6300,6364/[Ca II]447291,7324 line ratio of sev¬ 
eral core-collapse SNe of similar ages during nebular phase. 
Red symbols indicate SN type-11, orange type-llb, green type- 
lb, blue type-lc, and purple broad-lined type-lc. The SN phases 
with respect to the time of maximum light are shown next 
to the each data points and dashed vertical line indicates 
iPTF13bvn. The light-grey shaded region indicates the single 
star progenitors and the dark-grey shaded region indicates the 
binary progenitors, for type-lb/c SNe. The exact border be¬ 
tween these two regions is unknown, thus arbitrarily taken at 
[O l]/[Ca 11] = 1.5 in this plot for the purpose of indicating the 
two distinct p opulations. Refe rences fo r the spectra of individua l 


SNe: 19 87A (iPun et al 


1997D (iBenetti et al. 


19^: 1992H (IClocchiatti et al.i 


200 ih: 1998A (iPastorello et al.L 


1996); 


2005); 


120021 ) : 


1998bw IPatat et alX ^01 ): 1999em (Lemwde^all 
2004et ( Sahu et al. . 2006h ; 2004 gq (Maeda et all 2008h : 2005bf 
dMaeda et all l2007bl) : 2006br) (lOuimbv et al.L l2007h: 1996c b, 
2002ap, 2004fe, 2004gk, 2008ax. 2009if (iModiaz et all|20l4 . 


ISrivastav et al.L 120141) . In Figure^we show the early spectra of 
13bvn, around one week before and after maximum light com¬ 
pared to other SNe including 2007Y and 2010as. The absorption 
trough at 6200 A iniPTF13bv n has b een attributed to poss ibly 
Si 11 or Ne 1 bv ISrivastav et al.l (l2014l) and ICao et alJ (l2013h . In 
the case of SN 2010as, it was shown that this ab sorption is best 
matche d with high-velocity Ha rather than Si 11 (iFolatelli et al.L 
1201 4al) . 

3.1.4. H II region metallicity 

In the spectrum of iPTF13bvn there are narrow emission lines. 
The emission lines of Ha, [N ll]d6583, and [S 11]TT6716, 6731 
are presumably produced by the underlying H11 region at the SN 
position. Measurement of the logarithm of the flux ratio of [N11] 
to Ha emis sion lines, the N2 index , can be used as an indicator of 
metallicity. IPettini & Pagell (l2004l) gives the oxygen abundance 
as 12-i-log(0/H) = 8.90-1-0.57xA2.Our measurement of the A2 
index shows that the H 11 region associated with iPTF13bvn has 
the gas-phase metallicity of 12 -H log(0/H) = 8.63. This suggests 
that the metallicity of the progenitor o f iPTF13bvn is close t o the 
solar value of 12 -i- log(0/H) = 8.69 dAsnlund et al.L 1200^ . i.e. 
about 0.87 Z©. 

We also measured the oxygen abundance from a nearby H 
11 region 0.7 arcsec away from the SN, corresponding to pro¬ 
jected physical distance of 89 pc. We used the data taken as part 
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Fig. 8: Spectra of iPTF13bvn and SN 2008ax around Ha region. 
The flux of the SN 2008ax spectrum has been normalized and 
shifted. 


of the Carnegie Supernova Project (CSP; iHamuv et al.L l2006l) . 
using the 6.5 m Magellan Baade telescope at Las Campanas 
Obser vatory, Chile, and the instrument IMACS dBigelow et all 
Il998h on longslit spectroscopy mode with 300 lines/mm grating 
blazed at 4.3°. The observation was carried out in the night of 13 
August 2014 (local time), under 0.5 arcsec seeing. Measurement 
of the N2 index yields 12 H- log(0/H) = 8.62, which is very close 
to the derived on-site metallicity of 8.63. 

As already mentioned in subsection 13.1.21 at solar metal¬ 
licity the minimum initial mass for a single star to become a 
Wolf-Rayet star via stellar wind is ~25 M© dCrowtheii l2007h . 
Considering the metallicity estimate and core oxygen mass, it is 
unlikely that iPTF13bvn was produced by a >25 M© progenitor. 


3.2. Late-time light curve 


Our photometry unfortunately sampled only one epoch during 
the late-time evolution of iPTF13bvn. In Figure [10] we plot our 
UBVRc data points at 306 days together with the light curve 
constr ucted using published photometric data of ISrivastav et al.l 
(l2014t) . Assuming a constant decay rate since the latest points 
in the early light curves (-1-87 days in R and V, -1-52 days in 
B), the late-time tails of the light curves appear to decay more 
rapidly compared to t he ^®Co decay rate of 0.98 mag (100 d) ' 
(IWooslev et al.Lll989h . i.e. 1 .32 mag, 1.55 mag, and 1.13 mag 
(100 d)“' in R, V, and B bands respectiyely. 


The steeper light curye decay compared to ^®Co signi¬ 
fies that complete y-ray trapping did not occur in the case of 
iPTF13byn. This is a common phenomenon obseryed in other 
stripped-enyelope SNe dur ing the nebul ar phase e.g. SN 2008a x 
(lTaubenb erger_etaQ,|201T[). SN 1996N (ISollerman et aLLflOOSl) . 
SN 19901 dFlmhamdi et all l2004 . and is associated with rel- 
atiyely low ejecta mass. The decline of the light curye of SN 
19901 in nebular phase is the steepest in the V-band, followed 
by R and then B, which is likely also the case in iPTF13byn al¬ 
though the lack of data points during the late-time eyolution of 
the light curye preyents us from drawing such conclusion. 
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Rest —frame wavelength [A] 


Fig. 9; Early spectra of iPTFlSbvn at around one week before and after maximum light, compared to other stripped-envelope SNe. 
Vertical dashed lines show the positions of prominent helium lines, and dotted line sho w the position of high-vel ocity Ha. The 
spectra were taken from the Weizmann Interactive Supernova Data Repository (WISeREPi lYaron & Gal-YamLl2012ll 



Days after explosion 


Comparing the nebular [O I]/[Ca II] line ratio with other SNe, 
iPTF13bvn appears to be associated with low-mass progenitor 
stars. 

This piece of evidence provides another support to the sce¬ 
nario in which iPTF13bvn was produced by a hydrogen-poor 
progenitor star in a binary system that has undergone envelope 
stripping by means of close binary interaction. The initial mass 
of the progenitor star was not sufficiently high for the star to 
have become a Wolf-Rayet star, were it a single star. This re¬ 
sult shows very good agreement with the results from other 
methods in constraining the progenito r star, such as from hy- 
drody namical modeli n g JBersten et all l2014ll . analytical mod- 
eling dSrivastav et al.L 20141), and b i nary e volution modelling 
dBersten et al.L 2014 : lEldridge et all l2015h . In the context of 
hydrogen-poor SNe, iPTF13bvn further reinforces the impor¬ 
tance of massive close binaries as a prominent, if not the domi¬ 
nant, progenitors of these SNe. 


Fig. 10: Light curve of iPTF13bvn up until day 306 after the ex¬ 
plosion. Explosion time is i ndicated with dotted v ertical line. 
Early photometric data from ISrivastav et al.l (l2014l) are plotted 
in open squares; the latest points in the early light curves are 
connected to our data points (hlled squares) at nebular phase by 
dashed lines. Typical photometric error bars are smaller than the 
plotting symbol size. 


4. Summary 

We present late-time photometric and spectroscopic observa¬ 
tions of the type-Ib SN iPTF13bvn. A progenitor candidate of 
this SN has been detected in pre-explosion HST images, which 
is potentially the hrst for type Ib/c SNe. Our spectroscopy shows 
that the nebular spectrum of the SN is dominated by emis¬ 
sion lines of Mg I], [O I], and [Ca II]. Measurement of the 
flux of [O I]/l/16300, 6464 line shows that the oxygen mass of 
iPTF13bvn does not signihcantly exceed ~0.7 Mq, assuming 
that the emission is produced by 50% of the total oxygen mass, 
suggesting that the progenitor star could not have been more 
massive than ~15 - 17 Mq. The uncertainty originated from dis¬ 
tance estimate and oxygen line flux measurements are negligi¬ 
ble compared to the error from temperature and unseen oxygen. 
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